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We apply the continuous variable approach to study entangled dynamics of coupled 
harmonic oscillators interacting with a thermal reservoir and to a deterministic creation 
of entanglement in an atomic ensemble located inside a high-Q ring cavity. In the case 
of harmonic oscillators, we show that a suitable unitary transformation of the position 
and momentum operators transforms the system to a set of independent harmonic oscil- 
lators with only one of them coupled to the reservoir. Working in the Wigner represen- 
tation of the density operator, we find that the covariance matrix has a block diagonal 
form of smaller size matrices. This simple property allows to treat the problem to some 
extend analytically. We analyze the time evolution of an initial entanglement and find 
that the entanglement can persists for long times due to presence of constants of motion 
for the covariance matrix elements. In the case of an atomic ensemble located inside 
the cavity, the attention is focused on creation of one and two-mode continuous variable 
entangled states from the vacuum by applying laser pulses of a suitably adjusted ampli- 
tudes and phases. The pulses together with the cavity dissipation prepare the collective 
modes of the atomic ensemble in a desired entangled state. 
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1 Introduction 



Quantum information with continuous variables (CV) has attracted a great interest due 
to the simplicity in the generation, manipulation and detection of continuous variable states. 
Controlled dynamics and preservation of an initial entanglement encoded into CV states are 
challenging problems in quantum information technologies [ 1 1 and have led to the develop- 
ment of different techniques for generation, manipulation and detection of CV multipartite 
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entangled state [2|. It has been shown that CV entangled states have many applications 
in quantum information processing such as quantum teleportation 0, dense coding (ID, 
entanglement swapping [5 1, and quantum telecloning [6 1. 

The continuous variables approach is mostly applied to creation and manipulation of 
entangled states produced with squeezed light and linear optics. A wide theoretical atten- 
tion has been devoted to the application of the CV approach to the decoherence phenomena, 
in particular to the dissipative dynamics of open quantum systems, such as two coupled har- 
monic oscillators interacting with a dissipative reservoir. It has been analyzed with different 
approaches including the rotating-wave (RWA) and the Born-Markov approximations that 
assume a weak coupling between a system and the reservoir. A series of papers accounts 
these analysis for the case of two coupled harmonic oscillators interacting with a Marko- 
vian thermal reservoir and the work of Liu and Goan Q, Hu et al. fS], Maniscalco et 
al. (9) and Horhammer and Biittner ITOl accounts for a non-Markovian thermal bosonic 
reservoirs. Detailed discussions and extensive reference lists devoted to the decoherence of 
two harmonic oscillators can be found in Refs. [1 1-24]. In all these studies a general con- 
clusion made is that entanglement dynamics depends on the form of the reservoir and the 
non-Markovian nature of the reservoir preserves entanglement over a longer time. In this 
connection, we should mention work devoted to the variation of the entanglement sudden 
death phenomenon with the time scale of the evolution ||25H30ll . 

Recently, a good deal of attention has been given to the problem of application of 
the CV approach to creation and manipulation of entanglement in atomic systems |[3TH34ll . 
This is because a large collection of atoms can be efficiently coupled to quantum light and 
the existence of long atomic ground-state coherence lifetimes has been realized [35 1. Par- 
ticularly interesting are the studies of the dynamics of collective atomic spin states of driven 
two-level atoms and coupled to a cavity field [36]. It has been shown that the unconditional 
preparation of a two-mode squeezed state of effective bosonic modes can be realized in a 
pair of atomic ensembles interacting collectively with a two-mode optical cavity and laser 
fields [ 37 1 . Another studies considered the deterministic creation of cluster states between 
atomic ensembles ll38ll . 

In this paper, we illustrate applications of the CV approach to creation and processing 
of entanglement in two widely used systems: two coupled harmonic oscillators interacting 
with a non-Markovian thermal reservoir, and an ensemble of cold atoms located inside a a 
high-Q ring cavity. In the first case of the coupled harmonic oscillators, we introduce an 
unitary transformation of the position and momentum operators and find that in the trans- 
formed basis the system is represented by a set of independent harmonic oscillators with 
only one of them coupled to the environment. This fact makes the problem remarkably 
simple that the relaxation properties of coupled harmonic oscillators follow the same pat- 
tern as a single harmonic oscillator. We work within the correlation matrix representation, 
also known as the covariance matrix, and find that the transformation of the position and 
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momentum operators results in decoupled subsets of three and four equations of motion for 
the covariance matrix elements. This is particularly attractive for numerical analysis that it 
shorten the computations to the diagonalization or direct integration of small-size matrices. 
We then consider the time evolution of an initial entanglement encoded into the system of 
two coupled harmonic oscillators using the Gaussian continuous variable entangled states. 

In the second case, we consider a practical scheme for creation of entangled states 
in an ensemble of cold atoms located inside a a high-Q ring cavity. The scheme does 
not involve any external sources of squeezed light and networks of beam splitters used 
in the linear optics schemes l39l . It involves an atomic ensemble driven by laser pulses 
and CV entangled states are created by specifically chosen Rabi frequencies and phases 
of the pulses. The atoms interact with the laser pulses and the cavity field in a highly 
nonresonant dispersive manner that involves ground states, not the excited states of the 
atoms. Therefore, the process of creation of entangled states is not affected by the atomic 
spontaneous emission. 

2 Entangled dynamics of coupled harmonic oscillators 

We first consider entangled dynamics of a system composed of two mutually coupled 
identical harmonic oscillators of mass M and frequency 57 that are simultaneously inter- 
acting with a thermal reservoir. The dynamics are determined in terms of Gaussian con- 
tinuous variable entangled states, which are examples of multi-mode continuous variable 
entanglement states. We introduce a two-mode unitary transformation of the position and 
momentum operators and show that in the case of a non-Markovian reservoir, the formal- 
ism allows us to obtain, to some extend, an analytical solution for the dynamics of the 
system, and the exact analytical solutions under the Markov approximation. Thus, they are 
much more suited to give a united picture of this complex system. 

2.1 Hamiltonian of the system 

The system is determined by the Hamiltonian, which in terms of the position qi and 
momentum pi operators can be written as 

H = H S +H r + V ss +V sr , (2.1) 

where 

is the free Hamiltonian of the harmonic oscillators 

^=E(|t + ^") (2-3) 
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is the Hamiltonian of the reservoir, modeled as a set of independent oscillators of mass m n 
and frequency u„ to which the systems' oscillators are coupled, 

2 

v ss = \j2J2w (2 - 4) 

i=l j>i 

represents the interaction between the oscillators with the coupling strength A, and 

2 

Vsr = X X A «9n9i (2.5) 
n i—1 

represents a bilinear interaction between the oscillators and the reservoir with the 
strength A n . 

The Hamiltonian ( 12. It is written in terms of (qi,Pi) operators, usually called as the bare 
basis. However, this is not a very convenient basis and, as we shall see, the dynamics of 
the systems' harmonic oscillators is most transparently discussed in terms of transformed 
position and momentum operators that are obtained by an iV-mode unitary transformation 
of the position operators 02411401 

<7i = \J\ (qi - q_2) , <72 = (qi + 92) , (2.6) 

and the same transformation of the momentum operators 

Pi = (Pi -P2), P2 = (Pi + P2) ■ (2.7) 

The transformed operators satisfy the well known position-momentum commutation rela- 
tion, [q~i,pj] = ihSij. We note that the transformations from qi and pi to q~i and pi in- 
volve anti-symmetrical (qx,Pi) and symmetrical (q~2,P2) combinations of the position and 
the momentum operators, a close analog of the symmetric and antisymmetric multi-atom 
Dicke states B41II42I . 

In terms of the transformed operators, or equivalently in the transformed basis (q~i,pi), 
the Hamiltonian of the coupled oscillators H s + V s takes the form 



H s = H S + V S = J2(^7 + \Mtitql) , (2.8) 



where 



n 1 = n F = ^n 2 - A ^ = ^0? + -^, (2.9) 

and the interaction between the oscillators and the environment becomes 

V = V2j2^q n q2- (2-10) 
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One can see from Eqs. (12.8b and ( 12.10b that in terms of the transformed operators, the 
system is represented by two independent oscillators with only one being coupled to the 
environment. The oscillator effectively coupled to the environment is that one correspond- 
ing to the symmetric combination of the position and momentum operators. In addition, 
the effective frequency H2 of the oscillator coupled to the environment differs from f2i, the 
decoupled oscillator. Figure [2~T| illustrates the result of the unitary transformation of the 
system to new position and momentum operators. 

1X2 1 2 

Y\ ' 

(a) (b) 

Figure 2.1: A schematic diagram of a system of two coupled harmonic oscillators (a) before and (b) 
after the unitary transformation to new position and momentum basis. 

The oscillator effectively decoupled from the environment may be regarded as a 
decoherence-free oscillator. It should be stressed that the evolution of the oscillator is 
not free of the decoherence. We shall demonstrate that the oscillator still can evolve in time 
that may lead to decoherence. 




2.2 The master equation 

Let us now consider the evolution of the density operator of the system. We derive 
the master equation for the density operator and in the derivation we follow the standard 
method involving the Born approximation that corresponds to the second-order perturba- 
tive approach to the interaction between the oscillators and the environment, but we do 
not make the rotating-wave (RWA) and Markovian approximations. In the derivation, we 
assume that the system and the environment are uncorrected at t = 0, so that is the density 
operator of the total system factorizes pr(0) = p(0) (8> pr{0), where p(0) is the density 
operator of the harmonic oscillators, that we shall call the "system" oscillators, and p r (0) 
is density operator of the reservoir. Further, we consider the reservoir is in a thermal state 
of temperature T with the Boltzmann distribution of photons characterized by the mean 
occupation number 



N(u) = 




(2.11) 
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With these assumptions and after tracing over the reservoir modes, we arrive to the 
equation of motion for the reduced density matrix of the system of the form [7 8 1 

p(t) = - ~[# s + \Mnl{t)ql,p{t)] - ^ 2 (t)[q 2 ,{P2, P (t)}} 

- D 2 (t)[q 2 , [q 2 ,p(t)]} - ~/2(*)[g 2 , \p2,p(t)}]. (2.12) 

where represents a shift of the frequency of the oscillator due to the interaction with 
the environment, ^(t) is the dissipation coefficient, and D 2 (t) and f 2 (t) are diffusion 
coefficients. The parameters depend on the spectral density J(to) of the reservoir modes. 
They also depend on time, which results from the non-Markovian nature of the coupling of 
the oscillators to the reservoir. The explicit forms of the parameters are given in Refs. 08). 
As one could expect, the reservoir affects the evolution of only the oscillator 2 leaving the 
oscillator 1 to evolve freely in time. 

2.3 Covariance matrix of two harmonic oscillators 

In order to analyze the entangled dynamics of the oscillators, we solve the master equa- 
tion (I2.12l i. We apply the Wigner representation for the density operator. This requires a 
characteristic function, which can be written in terms of a covariance matrix, whose the 
elements are defined as 

V id = Tr | X - (AA^AX, + AXjAXi) p J , (2. 13) 

with AXi = Xi — (Xi) and X = q 2 ,p 2 , ■ ■ ■ , q~N iPn)- Using the definition ( 12.131 ) 

and the master equation (12.121 ) one can find the equations of motion for the covariance 
matrix elements that then can be solved for arbitrary initial conditions. The system of two 
harmonic oscillators is determined by a set of coupled linear differential equations for ten 
covariance matrix elements. In Refs. [7 8 1, the set of differential equations was solved using 
numerical methods. In what follows, we illustrate the advantage of working in the basis 
of the transformed position and momentum operators which will allow us to determine the 
covariance matrix elements in an effectively easy way requiring to solve separate sets of 
equations composed of three and four coupled differential equations. 

The set of inhomogeneous differential equations for the covariance matrix elements can 
be written in a matrix form as 

V(t) = A(t)V(t)+HF(t), (2.14) 
where we put the covariance matrix elements in a specific order in the column vector 



V(t) = co\{V 11 ,V 12 ,V 22 ,V 1 3,V u ,V 23 ,V 2 4,V 33 ,V 3 4,V44). (2.15) 
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This specific order of the elements Vij leads to a block diagonal form of the matrix of the 
coefficients A(t) with two blocks each involving three equations and one block involving 
four equations. The column vector of the inhomogeneous terms is of the form 

F(t) = co\(0,0,0,0,0,0,0,0,-f 2 (t),2hD 2 (t)). (2.16) 

The remaining elements are found from the symmetry property of the covariance ma- 
trix, Vij = Vji. 

From the set of three coupled equations of motion for the elements Vii, V± 2 , V 22 , one 
can find that there is a liner combination 

V+ = MQ%Vn + r^V 22 , (2.17) 

which equation of motion is Vq = and the remaining elements form a set of two coupled 
equations 

y-=40|Fi2, V n = -VCi, (2.18) 

where V{{ = MCl%Vn - {l/M)V 22 . 

The property of = indicates that the linear combination is a constant of 
motion, i.e. V{^(t) = V{^(0). In other words, V^(t) does not change in time and retains 
its initial value for all times. Physically, if initially the system was prepared in a state such 
that V^l (0) 7^ and with the other elements of the covariance matrix equal to zero, it would 
remain in that state for all times. For example, if the initial state is an entangled state, the 
initial entanglement of the system will remain constant in time. Therefore, the subspace 
composed of the V^[ (t) element can be regarded as a decoherence-free subspace. 

The remaining matrix elements V n and Vi 2 can undergo a temporal evolution. Since 
there is no damping involved in the equations of motion (12.181 , the solution would lead 
to the matrix elements continuously oscillating in time. It is easy to see, the solution of 
Eq. (12.181 1 has a simple form 

V n (t) = y n (0)cos(2ft F i) - 2O F y 12 (0)sin(2fi F t), 

V 12 (t) = V 12 (0)cos(2n F t) - ^P-sm(2n F t), (2.19) 

4s lp 

from which we see the matrix elements continuously oscillate in time with frequency 2flp. 
This indicates that the system will never reach a stationary time-independent state unless 
(0) = V12 (0) = 0. As we shall see, this feature will result in a continuous in time 
entanglement in the system. 

2.4 Initial two-mode squeezed state 

We have already seen that due to the presence of the constants of motion in the evolution 
of the covariance matrix elements, the dynamics of the system, even after a long time, may 
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strongly depend on the initial state. Since we are interested in the evolution of an initial 
two-mode entangled state and it is well known that two-mode squeezed states are examples 
of entangled states, we consider an initial squeezed vacuum state. We also demonstrate 
that with the specific initial state, the problem of treating the dynamics of two harmonic 
oscillators simplifies to analysis of the properties of only those constants of motion and the 
matrices which involve only the diagonal elements of the covariance matrix. As an initial 
state, we consider the following two-mode continuous variable squeezed state 

\ip )=e r ( h ^- b M)\O bl }\Q b2 ), (2.20) 

where r is the squeezing parameter and the kets 10^) and |0b 2 ) represent the state with zero 
photons in the modes 1 and 2, respectively. 

With the initial state (12.201 1. we easily find that non-zero initial covariance matrix ele- 
ments are 

V 1 i(0) = \ (Wc- 21 - ± ^e 2 ^ , V±(0) = ~ (Am%e 2r ± ^e" 2 ^ . (2.21) 

Note that the elements Vfi(0) and V^O) have opposite behavior with the squeezing cor- 
relations r. The frequency dependent part of V^(0) is squeezed while the other part is 
enhanced by the correlations. On the other hand, the frequency dependent part of V 3 ^ (0) is 
enhanced and the other part squeezed by the correlations. This is a crucial difference that 
will have a significant effect on the evolution of an entanglement injected into the system. 

2.5 Time evolution of an initial entanglement 

We now perform numerical analysis of time evolution of entanglement between two 
harmonic oscillators simultaneously coupled to a reservoir. We will illustrate the advan- 
tage of working in the transformed basis to obtain a simple interpretation of the results. 
In particular, to understand short time non-Markovian dynamics of entanglement and to 
provide conditions for optimal and stable long time entanglement. We adopt the negative 
partial transpose criterion that is known as the necessary and sufficient condition for en- 
tanglement of a two-mode Gaussian state [45 , 46 1. We will use this criterion to quantify 
the amount of entanglement in the system, and will denote it by a parameter . In all 
examples considered here, we assume that the oscillators interact with a reservoir of tem- 
perature T = lOhn/ks and were initially in the two-mode squeezed state \ipo). Moreover, 
we assume a Gaussian-type spectral density for the reservoir modes 

^^^.V/A. (2 . 22) 

where M is the mass of system, A is cut-off frequency that represents the highest frequency 
in the reservoir, 70 is proportional to the coupling strength between the oscillators and the 
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Figure 2.2: Time evolution of the negativity r\- for 70 = 0.1, A = 100, A = and different r: 
r = 1.0 (solid line), r = 1.498 (dashed line), r = 2.0 (dashed-dotted line). 



environment and n determines the type of the reservoir. We choose n = 1, the so-called 
Ohmic reservoir. 

We first consider the case of mutually independent oscillators with A = 0, but coupled 
to a thermal reservoir. The interaction with the reservoir is usually recognized as an irre- 
versible loss of entanglement and information encoded in the internal states of the system 
and thus is regarded as the main obstacle in practical implementation of entanglement. Fig- 
ure |22] shows the negativity as a function of time for the initial two-mode squeezed state 
with different degree of squeezing r. First of all, we note that there is a threshold value 
for the degree of squeezing r at which a continuous in time entanglement occurs 1 19 1. The 
threshold that corresponds to entanglement undergoing the phenomenon of sudden death, 
occurs at r = 1.498. 

The presence of the threshold value for r at which continuous in time entanglement 
occurs has a simple interpretation in terms of the covariant matrix elements. Consider the 
threshold in the long time limit in which we may consider the evolution under the Markov 
approximation, but retaining the non-RWA terms. Under this approximation, we can put 
72 (t) — > 70 which then allows us to obtain a simple analytical solution for the threshold 
condition for entanglement. It is easy to show that the threshold for two mode entanglement 
occurs at 

V11 (t)V«(t) = 1, (2.23) 

so that the two modes are entangled when Vn {t)Vi4,(t) < 1, otherwise are separable. 

Note that the covariance matrix element Vn(t) is associated with the relaxation free 
modes whereas the element V44 (t) is associated with the mode that is coupled to the reser- 
voir and thus undergoes the damping process. Under the Markov approximation, we find 
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from Eq. (12.14b that in the long time limit of t 3> 7 2 , the element Va{t) reaches the 
stationary value equal to the level of the thermal fluctuations 

Vu(t)^-2N+1, (2.24) 

whereas Vn(t) retains its time dependent behavior which depends on the initial values 

V U ® = V u (0) cos 2 n F t + ^ (^f) 2 1 2 - (2-25) 

We point out that the dependence on the initial values of the long time behavior of Vn(t) 
is due to the presence of the constant of motion Vyy. 

Averaging Eq. ( |2.25t over a long period of oscillations, the threshold condition ( 12.231 ) 
simplifies to 

2Vn(0) (2JV + 1) =1. (2.26) 

We see that the threshold behavior of entanglement depends on the initial value of the 
covariance matrix element Vn (0). In other words, the entanglement behavior can be con- 
trolled by the suitable choosing of the initial state. For example, with the initial state ( |2.20| i, 
Vii (0) = exp(— 2r) /2, and then we find from Eq. (12. 26b that continuous entanglement oc- 
curs for the degree of squeezing 

r > - In (2N + l) . (2.27) 

With the parameter value T = lOMl/fcs, we find that the threshold value for r equals 
to 1.498 that is the same found numerically and plot in Fig. 12.21 We should point out here 
that the same condition for the threshold value of r has been found under the RWA approx- 
imation fl9l . Thus, we may conclude that the threshold value for continuous entanglement 
is not sensitive to the RWA approximation. 

The preservation of an initial entanglement over a long time is usually related to the 
non-Markovian nature of the reservoir. We stress that the continuous in time oscillations 
are related to the Markovian rather than to the non-Markovian nature of the reservoir as the 
matrix elements Vii (t) and V22 (t) determine dynamics of the oscillator that is not coupled 
to the reservoir. To show this more qualitatively, we plot in Fig. [23] the time evolution of 
an initial entanglement for different A, corresponding to the bandwith of the reservoir. It is 
evident from the figure that the amplitude of the oscillations increases with increasing A in- 
dicating that the oscillations are associated with broadband rather than narrow-band nature 
of the reservoir. 

It is also interesting to discuss the dependence of the long time entanglement on the 
relaxation rate 70. An example of this feature is shown in Fig. 12.41 It is interesting to note 
that under the relaxation the entanglement oscillates in time and the amplitude of the oscil- 
lations increases with increasing 70 leading to a better entanglement when the oscillators 
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Figure 2.3: Time evolution of the negativity r/_ for A = 0, r = 1.6, 70 = 1.0 and different values of 
the cut-off frequency A :A = 200 (solid line), A = 500 (dashed line), A = 800 (dashed-dotted line). 

are strongly damped. It is a surprising result as one could expect that entanglement should 
decrease with increasing 70. Again, a straightforward interpretation of this effect can be 
gained from a qualitative inspection of the properties of the transformed covariance matrix. 




Figure 2.4: Time evolution of the negativity r\- for A = 100, A = 0, r = 1.6 and different values of 
the relaxation rate 70: 70 = 0.05 (solid line), 70 = 1.0 (dashed line), 70 = 5.0 (dashed-dotted line). 

It is easily verified that in the limit of vanishing damping, 72 (t) — > and A = 0, 
the equations of motion of the matrix elements V33 , V34 , V44 reduce to that of the matrix 
elements Vu, V12, Vii- One could argue that in this limit the covariance matrix elements 
V33, V34, V44 coincidence with the elements Vu, V12, Vii- Of course, their time evolution 
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is determined by the same equations, but there is a subtle difference in their initial values. 
For example, the initial values of the linear combinations VfJ are 

ttf(O) = \ (W F e" 2r ± ^e 2 ^ , (2.28) 

whereas 

V£(P) = \ (M£lW r ± ^- 2r ) ■ (2.29) 

The initial elements are significantly different that what appears as a squeezed component 
in Vff(O), the counterpart in V^O) appears as an anti-squeezed component. This is a cru- 
cial difference that has a significant effect on the evolution of an entanglement. These two 
contributions cancel each other that results in no oscillations in the entanglement evolution 
when 70 <C 1. On the other hand, for large 70 the covariance matrix elements V33, V34, V44 
are rapidly damped to their stationary values leaving the elements V\\, V12, V22 continu- 
ously oscillating in time. These oscillations lead to continuous oscillations of the entangle- 
ment seen in Fig. 12.41 

The same arguments apply for the presence of the oscillations in the limit of a large A, 
seen in Fig. 12.31 Using similar arguments, one can finds that for large A, the covariance 
matrix elements V33 , V34 , V44 are rapidly damped to their stationary values leaving the 
elements Vn, V12, V22 continuously oscillating in time. 
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Figure 2.5: Time evolution of the negativity r\- for 70 = 0.1, A = 100, A = 0.8 and different r: 
r = 1.0 (solid line), r = 1.498 (dashed line), r = 2.0 (dashed-dotted line). 

Figure|23]shows the evolution of entanglement when the oscillators are coupled to each 
other. In this case there is no continuous stationary entanglement. Thus, the interaction 
between the oscillators has a destructive effect on the stationary entanglement. However, 
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for a large squeezing, entanglement re-appears in some discrete periods of time, exhibiting 
periodic sudden death and revival of entanglement. In other words, the threshold behavior 
of entanglement is a periodic function of time. 

As before, this feature has a simple interpretation in terms of the covariance matrix el- 
ements. According to Eq. ( I2.23K for a given temperature the threshold value for entangle- 
ment depends on the covariance matrix element Vu(t) which, on the other hand, depends 
on A through the frequency parameter Hp- We see from Eq. ( 12.91 ) that the frequency ftp 
decreases with increasing A. Thus, according to Eq. (12.25b for interacting oscillators the 
matrix element Vu(t) oscillates slowly in time. In this case, the averaging over the oscil- 
lations is not justified and thus the threshold condition for entanglement is the oscillating 
function of time even in a long time regime. 

3 Creation of CV entangled states in an atomic ensemble 

In this section, we consider a completely different system that addresses the issue of the 
creation of CV entangled states in an atomic system. The system we consider is composed 
of an atomic ensemble located inside a high-Q ring cavity with two mutually counter- 
propagating modes, as illustrated in Fig. 13.61 In this two-mode configuration, the photon 




Figure 3.6: A schematic diagram of a two-mode ring cavity containing an ensemble of cold atoms 
trapped along the cavity axis. The cavity modes are damped with the same rate k. The driving laser 
fields of the Rabi frequencies Q u , Q a and phases (j> u , <j> s are injected through the cavity mirrors and 
co-propagate with one of the cavity modes. 

redistribution can occur between the counter-propagating modes which, as we shall see, 
may result in multi-mode squeezing. The procedure is associated with the collective dy- 
namics of the atomic ensembles subjected to driving lasers of a suitably adjusted amplitudes 
and phases. 
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Recently, Parkins et al. 071 have proposed a scheme for preparing single and two-mode 
squeezed states in atomic ensembles located inside a single mode high-Q ring cavity. The 
scheme is based on a suitable driving of the atomic ensembles with two external laser fields 
that prepares the atoms in a pure squeezed (entangled) state. Similar schemes have been 
proposed to realize an effective Dicke model operating in the phase transition regime |47| . 
to create a stationary subradiant state in an ultracold atomic gas ||48]]. This approach has 
recently been proposed as a practical scheme to prepare trapped and cooled ions in pure 
entangled vibrational states [49 1 and to prepare four ensembles of hot atoms in pure entan- 
gled cluster states [38). The scheme proposed by Parkins et al. 11371 assumes that atomic 
ensembles are effectively coupled to a single or two modes co-propagating with the driving 
fields. In fact, a practical ring cavity is composed of two degenerate mutually counter- 
propagating modes that can simultaneously couple to the atomic ensembles with the same 
coupling strengths [50- 52]. Here, we explore the role of the counter propagating mode in 
the creation of one and two-mode entangled states in a single atomic ensemble. 



3.1 Hamiltonian of the system 

We consider a system composed of a single atomic ensemble located along the axis of 
a high-Q ring cavity, as illustrated in Fig. 13.61 The cavity is composed of three mirrors 
that create two mutually counter-propagating modes, called clockwise and anti-clockwise 
modes, to which the atoms are equally coupled. The cavity modes are degenerate in fre- 
quency, i.e. oj + = tjj- = oj c and are damped with the same rate k. The atomic ensemble 
contains a large number of identical four-level atoms, each composed of two ground states, 
\0j), \lj), and two excited states \uj), \sj), where the subscript j labels atoms in the en- 
semble. The ground state |0 3 -) of energy Eq = is coupled to the excited state \sj) by 
a laser field of the Rabi frequency £l s and frequency ojl s that is detuned from the atomic 
transition frequency by A s = (lj s — Wi s ), where u s = E s /h and E s is the energy of 
the state \sj). Similarly, the ground state of energy E\ = f\W\ is coupled to the ex- 
cited state \uj) of energy E u — ftw u by an another laser field with the Rabi frequency Cl u , 
and the angular frequency ojl u that is detuned from the atomic transition — > \uj) by 
A u = (lu u — lo\ — LdLu)- The frequencies u>Lu and ujls of the laser fields are matched 
close to the cavity frequency, so that the wave numbers of the laser fields k u and k s are 
approximated by k u ss k s sa k. 

The Hamiltonian of the system, in the rotating-wave approximation and in the interac- 
tion picture, can be written as 



H — H + Hal + Hao 



(3.1) 
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where 

N 

H = »5Z(w u |uj)<Uj|+w a |sj)<«j|+wi|lj)<lj|) 

3=1 

+ tvjj c {a [ + a + + ala_) (3.2) 
is the free Hamiltonian of the atoms and the cavity modes, 

1 N 

Hal = -h^2{si u e- i ^ t+ ^- k ^\u j )(l j \ 

3 = 1 

+ n a e- i [( wl '* + *^- fc '^]|s J -)<0 J -|+H.c.} (3.3) 
is the interaction Hamiltonian between the atoms and the driving laser fields, and 

N 

H AC = fi£ {ffi +) e^ K)(0>+ + gi-h- ikx > h)<0>_ 

+ 5 ( + )e^| S ,)(l j |a ++5 (-)e-^> j )(l j |a_+H.c.} (3.4) 

is the interaction Hamiltonian between the atoms and the cavity modes. The notation used 
here should be interpreted as follows. The operators a± and aj_ are the annihilation and 
creation operators of the two counter-propagating modes of the cavity; clockwise and anti- 
clockwise propagating modes, and k is their wave number. The parameters and gi^ 
are the coupling constants of the atomic transitions to the cavity modes, fi u and Q s are 
Rabi frequencies of the driving laser fields, and (fi u , <p s are their phases. We have set the 
energy of the ground state |0j) equal to zero and have denoted the energies of the excited 
atomic levels by fkdi (i = 1, u, s). 

We now introduce detunings of the laser fields from the atomic transition frequencies 

A„ = (uj u - wi) - ui Lu , A s =lj s -uj Ls , (3.5) 

and assume that the laser frequencies satisfy the resonance condition uil s — ujl u = 2u!\. 
Next, we make few standard approximations on the Hamiltonian to eliminate the excited 
states to neglect atomic spontaneous emission and to obtain an effective two-level Raman- 
coupled Hamiltonian. In order to eliminate spontaneous scattering of photons to modes 
other than the privileged cavity mode, we assume that the detunings are much larger 
than the Rabi frequencies, cavity coupling constants and the atomic spontaneous emission 
rates, i.e. 

A„, A s > gW, g&\ n a , ft u , 7a , 7u, (3.6) 

where 7 S and j u are the total spontaneous emission rates from the states \sj) and \uj), 
respectively. The assumption about large detunings allows us to perform the standard adi- 
abatic elimination of the atomic excited states, \sj), \uj), and obtain an effective two-level 
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Hamiltonian involving only the ground states of the atoms 





' 9l 




( 9l _ 


.9? 













(a^a + + a^a_^j 



-4= { Pu (j *4 + J 2 feaL) + Pa (4*°+ + J -2fc a -) + H - c -} . ( 3 - 7 > 



iV 

where we have assumed that the coupling constants g u = and g s = <?g are the same 
for both modes, S c = uj c — {u>l s — cj%) is the detuning of the cavity frequency from the 
Raman coupling resonance, 

N N 

J* = 2 E (l^)^l - > Jmk = J2 Wile*"*"* (3-8) 

are position dependent collective atomic operators, and 

A.=VjV^e-^, /3 s = ViV^e-^ (3.9) 

are the coupling strengths of the effective two-level system to the cavity modes. 

The first line in Eq. ( 13.7b represents the free energy of the atomic ensemble and the 
intensity dependent (Stark) shift of the atomic energy levels. The second line represent the 
interaction between the cavity field and the atomic ensemble. The three collective atomic 
operators, Jofc, Jik and J -2k, which appear in Eq. (13.71 i. arise naturally for the position 
dependent atomic transition operators and appear in a cavity with two mutually counter- 
propagating modes. In the case of a single-mode cavity, the Hamiltonian involves only 
the J ofc operator [|37ll47ll48l . 

To avoid unessential complexity due to the presence of the free energy and the Stark 
shift terms in the Hamiltonian (13.71 >. we will work with a simplified version of the Hamil- 
tonian by choosing the frequencies of the driving and the cavity fields such that 

- i-i,- 

With this choice of the parameters, the Hamiltonian ( 13. 7t reduces to 

H e = -j= { P u (Jofcfl^ + Jafcffll) + Pa (4*4 + J-2k a -) + H ' C -} ■ ( 3 - U ) 



h 

N 

The parameters of the Hamiltonian are a function of the detunings and Rabi frequencies 
of the two highly detuned laser fields co-propagating with the clockwise cavity mode and 
thus could be controlled through the laser frequencies and intensities. 

Finally, we reformulate the Hamiltonian J3.1 It in terms of bosonic variables by adopt- 
ing the Holstein-Primakoff representation of angular momentum operators ll53ll . In this 
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representation, the collective atomic operators, J mk , J m k and J z are expressed in terms of 
annihilation and creation operators C m k an d Clfc °f a sm gl e bosonic mode 



J: 



ink 



C mk y/ N - Cl k C mk , J z =Cl k C mk -N/2. (3.12) 



Provided the atoms in each ensemble are initially prepared in their ground states 1 j ) , and 
taking into account that due to large detunings of the driving fields, the excitation proba- 
bility of each atom is low during the laser-atom-cavity coupling, i.e., the number of atoms 
transferred to the states \lj) is expected to be much smaller than the total number of atoms 
in each ensemble, i.e., (Cj nk C mk ) <C N. By expanding the square root in Eq. (I3.12l i and 
neglecting terms of the order of 0(1/ N), the collective atomic operators can be approxi- 
mated as 



Jmk = VNC mk , J z = -N/2, (3.13) 



where 



C mk = -j= V b je imkx * , to = 0, ±2, (3.14) 

are collective bosonic operators with the operators bj and fcj obeying the standard bosonic 
commutation relation = 5 ji- 

lt is easy to prove that in the limit of N S> l,we obtain the standard bosonic commuta- 
tion relation 



(3.15) 



Thus, in terms of the collective bosonic operators Co k and C±2 k , the effective Hamil- 
tonian ( 13. lit takes the form 



h (f3uC ok + p s Cl k ) al+h (PuC 2k + PsCl 2k ) al + H.c. (3.16) 



The important property of the bosonic representation is the fact that the Hamiltonian of an 
ensemble of atoms trapped inside a ring cavity can be expressed as the interaction between 
the cavity modes and three orthogonal field modes. 

Since our objective is to prepare the atomic ensemble in a desired entangled state in the 
presence of a possible loss of photons due to the damping of the cavity mode, we shall work 
with the density operator p of the system whose the evolution is governed by the master 
equation 



where 



P =- l -[H e , P ]+C cP , (3.17) 
n 

C c p = 2 K ^2 ( 2ai /°°I _ a l a iP - P a l a i) j (3.18) 
»=± 
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is the Liuivilian operator representing the damping of the cavity field modes with the rate k. 

In what follows, we show that by a proper choosing of the Rabi frequencies fi u , H s of 
the laser fields and their phases, (j> u , <f> s one may prepare the system in a desired state that 
then decays with the rate n to a steady-state form. 



3.2 One and two-mode entangled states 

Let us now illustrate in details how to construct entangled states in a single ensemble of 
cold atoms located inside a two-mode ring cavity. We focus on the creation of single and 
two-mode entangled states and attempt to characterize the entanglement in terms of unitary 
operators known as squeezed operators. The operators for single and two-mode squeezed 
states are defined as B54II551 



So (Co) = exp 



2 {^° C o" k ~ 



Sfcfcfr) = exp fec 2 fcC_ 2fc - CxC^cO , (3.19) 



2fc v - / -2fc J ' 

where £o an d £i are complex one and two-mode squeezing parameters, respectively. 

Entangled (squeezed) states in a single ensemble are constructed from the vacuum by 
a unitary transformation associated with the realistic dynamical process determined by the 
master equation ( 13.16b . The creation of entangled states is done in two steps. In the first 
step, we adjust the driving lasers to propagate in the clockwise direction, along the cav- 
ity mode (+). In this case, the dynamics of the system are determined by the Hamilto- 
nian ( 13.161 ). We then send series of laser pulses of phases <j> u — <f) s — and arbitrary Rabi 
frequencies fl u and f2 s , but such that \f3 u \ > \/3 s \. With this choice of the parameters of the 
driving lasers and under the unitary squeezing transformation 

So(-£o)S±fc(-£i)pSo(£o)S±fc(£i) = P, (3.20) 

with 

the master equation ( 13.171 ) becomes 

j t P = -i[H e ,p\+C c p, (3.22) 



where 



H e - So(-£o)S±fc(-£i)ir e So(£o)S± fc (£i) 

= M/3„| 2 - W (a\C 0k + alC 2k + H.c.) . (3.23) 
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It is seen that under the squeezing transformation, the Hamiltonian represents a simple 
system of two independent linear mixers, where the collective bosonic modes Cofe and C2k 
linearly couple to the cavity modes a + and a_, respectively. The mode C-2fc is decoupled 
from the cavity modes and therefore does not evolve. In other words, the state of the 
mode C-2k cannot be determined by the evolution operator for the Hamiltonian ( 13.231 ). 
The important property of the transformed system is that the master equation ( 13.221 ) is 
fully soluble, i.e. all eigenvectors and eigenvalues can be obtained exactly. Hence, we can 
monitor the evolution of the bosonic modes towards their steady-state values. Since we are 
interested in the steady-state of the system, we confine our attention only to the eigenvalues 
of Eq. ( 13.22b . which are of the form 

(3.24) 

Evidently, both eigenvalues have negative real parts which means that the system subjected 
to a series of laser pulses up to a short time t will evolve (decay) to a stationary state that 
is a vacuum state. Thus, as a result of the interaction given by the Hamiltonian ( 13.231 ), and 
after a sufficiently long evolution time, the modes a±, Cofc and C2k will be found in the 
vacuum state, whereas the mode C-2k will remain in an undetermined state. The state of 
the mode C_2fc will be determined in the next, second step of the preparation process. 

In order to estimate the time scale for the system to reach the steady-state, we see 
from Eq. ( 13.24b that as long as ^/|/3„| 2 — \(3 S \ 2 > k/2, the time scale for the system to 
reach the steady state is of order of ~ 2/k. Thus, as a result of the cavity damping the 
system, after a sufficient long time, will definitely be found in the stationary state. 

Thus, after the first step of the preparation, we find that in the steady-state the density 
matrix representing the state of the transformed system is in the factorized form 

p(r ~ 2/k) = p v ® p C -2k ) (3-25) 

where 

Pv = |0 a+ ,0 a _,0c Ofc ,0 C2fc )(0 a+ ,0 a _,0 Cofc ,0c 2 J (3.26) 

is the density matrix of the four modes prepared in their vacuum states, and pc- 2 k * s me 
density matrix of the mode C_2fc whose the state has not yet been determined. The ket 
1 0a + 1 a _ , 0c Ofe , 0c 2fc ) represents the state with zero photons in each of the modes. 

What left is to prepare the remaining collective mode C_2fc in a desired squeezed vac- 
uum state. This is done in what we call the second step of the preparation, in which we 
first adjust the driving lasers to propagate along the anti-clockwise mode (— ). We then 
send series of pulses of frequencies, phases and amplitudes the same as in the above first 
step. As a result of the coupling to the cavity mode (— ), the interaction is now governed 



»?± = -2 ± 
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by the Hamiltonian ( 13.18b . and therefore after the unitary squeezing transformation the 
Hamiltonian of the system takes the form 

He = M-£o)S ±fe (-£i)ff e S (£o)S±fe(&) 

= VKI 2 - W ( a - C ™ + alC- 2k + H.c.) . (3.27) 

As above in the case of the coupling to the cavity mode a + , the Hamiltonian d3.27| i 
describes a system of two independent linear mixers. Hence, the state of the system 
will evolve during the interaction towards its stationary value, and after a suitably long 
time, ~ 2/k, the transformed system will be found in the vacuum state. 

Hence, after the second step of the preparation, the transformed system is found in the 
pure vacuum state determined by the density matrix of the form 

p(r~ 4//s) = |*><*l> (3.28) 

where 

I*) - So{-to)S± k (-Zi)\9) = |0a + ,0 a _,0c Ofe ,0c 2fc ,0 c _ 2 J (3.29) 

represents the vacuum state of the transformed system and the ket \^) represents the final 
stationary state of the system. 

If we now perform the inverse transformation from \ty) to |<3>), we find that the system 
is in the multi-mode pure squeezed state 

I*) = S (Zo)S±k(Zi)\ac oh ,Oc 2k ,Oc_ 2k ) ® |0 a+ ,0 Q _). (3.30) 
The density operator representing the multi-mode squeezed state ( 13.301 ) is 

P = So(£o)S± fc (£i)|{0})({0}|So(-£o)S±fc(-a), (3.31) 

where |{0}) = |0 Cofc , 0c 2k , G _ 2h ) ® |0 a+ ,0 Q _). 

Equations ( 13.30b and (I3.31l l show that in the steady-state the cavity modes are left in 
the vacuum state and the atomic ensemble is prepared in single and two-mode squeezed 
states. In other words, it shows that the collective mode Cofc is prepared in the one- 
mode squeezed vacuum state Sb(£o)|0cofc)' whereas the collective modes C±2k are in the 
two-mode squeezed vacuum state S±k{£,i)\0c 2k , 0c_ 2fc ) associated with the superposition 
of two, position dependent counter-propagating modes. The superposition nature of the 
squeezed states allows us to conclude that the atomic ensemble, after the interaction with 
the sequences of the laser pulses is prepared in single and two-mode entangled states. 

4 Conclusions 

We have considered to examples of application of the CV approach to entanglement 
creation and processing. In the first, we have analyzed the time evolution of an initial 
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entanglement encoded into two coupled harmonic oscillators interacting with a thermal 
reservoir. We have shown that it is very convenient to analyze dynamics of the oscillators 
in terms of the evolution of the covariance matrix elements. By performing a suitable 
transformation of the position and momentum operators of the system oscillators, we have 
shown that the set of ten coupled differential equations for the covariance matrix elements 
splits into decoupled subsets of smaller sizes involving only three and four equations. A 
general feature of the entanglement evolution is that it exhibits two characteristic time 
scales, a shot time regime where an initial entanglement is rapidly damped and a long time 
regime where the entanglement undergoes continuous undamped oscillations. Depending 
on the initial amount of entanglement encoded into the system, it can be preserved for all 
times or may undergo the sudden death and revival phenomena. We have also found that 
in contrast to what one could expect, a stronger coupling of the oscillators to the thermal 
reservoir leads to a better stationary entanglement than in the case of a weak coupling. 

In the second example of application of the CV approach, we have described a practical 
scheme for the creation of CV entangled states of effective bosonic modes realized in an 
atomic ensemble interacting collectively with two counter propagating modes of a ring 
cavity. The basic idea of the scheme is to transfer the ensemble field modes into suitable 
linear combinations that can be prepared, by a sequential application of the laser pulses, in 
pure squeezed vacuum (entangled) states. 

References 

[1] S. L. Braunstein and A. K. Pati, Quantum Information Theory with Continuous Vari- 
ables, Kluwer, Dordrecht, 2003. 

[2] S. L. Braunstein and P. van Loock, Quantum information with continuous variables, 
Rev. Mod. Phys. 77 (2005), 513-577. 

[3] A. Furusawa, J. L. Sorensen, and S. L. Braunstein, C. A. Fuchs, and H. J. Kimble and 
E. S. Polzik, Unconditional quantum teleportation, Science 282 (1998), 706-709; N. 
Takei, T. Aoki, S. Koike, K. Yoshino, K. Wakui, H. Yonezawa, T. Hiraoka, J. Mizuno, 
M. Takeoka, M. Ban, and A. Furusawa, Experimental demonstration of quantum tele- 
portation of a squeezed state, Phys. Rev. A 72 (2005), 042304; H. Yonezawa, S. L. 
Braunstein, and A. Furusawa, Experimental demonstration of quantum teleportation 
of broadband squeezing, Phys. Rev. Lett. 99 (2007), 1 10503. 

[4] X. Y. Li, Q. Pan, J. T. Jing, J. Zhang, C. D. Xie, and K. Peng, Quantum dense cod- 
ing exploiting a bright Einstein-Podolsky-Rosen beam, Phys. Rev. Lett. 88 (2002), 
047904. 

[5] X. J. Jia, X. L. Su, Q. Pan, J. R. Gao, C. D. Xie, and K. Peng, Experimental demon- 
stration of unconditional entanglement swapping for continuous variables, Phys. Rev. 
Lett. 93 (2004), 250503. 



22 



Li-hui Sun et al. 



[6] N. Takei, H. Yonezawa, T. Aoki, and A. Furusawa, High-fidelity teleportation beyond 
the no-cloning limit and entanglement swapping for continuous variables, Phys. Rev. 
Lett. 94 (2005), 220502; J. Zhang, C. Xie, K. Peng, and P. van Loock, Anyon statistics 
with continuous variables, Phys. Rev. A 78 (2008), 052121. 

[7] K.-L. Liu and H.-S. Goan, Non-Markovian entanglement dynamics of quantum con- 
tinuous variable systems in thermal environments, Phys. Rev. A 76 (2007), 022312. 

[8] B. L. Hu, J. P. Paz and Y. Zhang, Quantum Browian-motion in a general environment 
- exact master equation with nonlocal dissipation and colored noise, Phys. Rev. D 45 
(1992), 2843-2861. 

[9] S. Maniscalco, S. Olivares, and M. G. A. Paris, Entanglement oscillations in non- 
Markovian quantum channels, Phys. Rev. A 75 (2007), 0621 19. 

[10] C. Horhammer and H. Biittner, Environment-induced two-mode entanglement in 
quantum Brownian motion, Phys. Rev. A 11 (2008), 042305. 

[11] D. Braun, Creation of entanglement by interaction with a common heat bath, Phys. 
Rev. Lett. 89 (2002), 277901. 

[12] M. S. Kim, J. Lee, D. Ahn, and P. L. Knight, Entanglement induced by a single-mode 
heat environment, Phys. Rev. A 65 (2002), 040101. 

[13] F. Benatti, R. Floreanini, and M. Piani, Environment induced entanglement in Marko- 
vian dissipative dynamics, Phys. Rev. Lett. 91 (2003), 70402. 

[14] S. Oh and J. Kim, Entanglement between qubits induced by a common environment 
with a gap, Phys. Rev. A 73 (2006), 062306. 

[15] B. Bellomo, R. Lo Franco, and G. Compagno, Non-Markovian effects on the dynam- 
ics of entanglement, Phys. Rev. Lett. 99 (2007), 160502; Entanglement dynamics of 
two independent qubits in environments with and without memory, Phys. Rev. A 11 
(2008), 032342. 

[16] M. G. A. Paris, Optimized quantum nondemolition measurement of a field quadrature, 

Phys. Rev. A 65 (2002), 012110. 
[17] A. Serafini, F. Illuminati, M. G. A. Paris, and S. De Siena, Entanglement and purity 

of two-mode Gaussian states in noisy channels, Phys. Rev. A 69 (2004), 022318. 
[18] P. J. Dodd and J. J. Halliwell, Disentanglement and decoherence by open system 

dynamics, Phys. Rev. A 69 (2004), 052105. 
[19] J. S. Prauzner-Bechcicki, Two-mode squeezed vacuum state coupled to the common 

thermal reservoir, J. Phys. A: Math. Gen. 37 (2004), L173-L181. 
[20] F. Benatti and R. Floreanini, Entangling oscillators through environment noise, J. 

Phys. A: Math. Gen. 39 (2006), 2689-2699. 
[21] J.-H. An and W.-M. Zhang, Non-Markovian entanglement dynamics of noisy 

continuous-variable quantum channels, Phys. Rev. A 76 (2007), 042127. 
[22] C.-H. Chou, T. Yu, and B. L. Hu, Exact master equation and quantum decoherence of 

two coupled harmonic oscillators in a general environment, Phys. Rev. E 11 (2008), 

011112. 



Continuous variables approach to entanglement 



23 



[23] M. Paternostro, G. Adesso, and S. Campbell, Passing quantum correlations to qubits 
using any two-mode state, Phys. Rev. A 80 (2009), 062318. 

[24] J. P. Paz and A. J. Roncaglia, Dynamics of the entanglement between two oscillators 
in the same environment, Phys. Rev. Lett. 100 (2008), 220401; Dynamical phases 
for the evolution of the entanglement between two oscillators coupled to the same 
environment, Phys. Rev. A 79 (2009), 032102. 

[25] T. Yu and J. H. Eberly, Finite-time disentanglement via spontaneous emission, Phys. 
Rev. Lett. 93 (2004), 140404. 

[26] T. Yu and J. H. Eberly, Sudden Death of Entanglement, Science 323 (2009), 598-601. 

[27] Z. Ficek and R. Tanas, Dark periods and revivals of entanglement in a two-qubit 
system, Phys. Rev. A 74 (2006), 024304. 

[28] X. Cao and H. Zheng, Non-Markovian disentanglement dynamics of a two-qubit sys- 
tem, Phys. Rev. A 77 (2008), 022320. 

[29] L. Mazzola, S. Maniscalco, J. Piilo, K.-A. Suominen, and B. M. Garraway, Sudden 
death and sudden birth of entanglement in common structured reservoirs, Phys. Rev. 
A 79 (2009), 042302. 

[30] J. Jing, Z.-G. Lii, and Z. Ficek, Breakdown of the rotating-wave approximation in 
the description of entanglement of spin-anticorrelated states, Phys. Rev. A 79 (2009), 
044305. 

[31] M. D. Lukin, Colloquium: Trapping and manipulating photon states in atomic ensem- 
bles, Rev. Mod. Phys. 75 (2003), 457-472. 

[32] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Electromagnetically induced 
transparency: Optics in coherent media, Rev. Mod. Phys. 77 (2005), 663-673. 

[33] K. Hammerer, A. S. Sorensen, and E. S. Polzik, Quantum interface between light and 
atomic ensembles, Rev. Mod. Phys. 82 (2010), 1041-1093. 

[34] G.-X. Li and Z. Ficek, Creation of pure multi-mode entangled states in a ring cavity, 
Optics Commun. 283 (2010), 814-821. 

[35] D. N. Matsukevich, T. Chaneliere, S. D. Jenkins, S. Y. Lan, T. A. B. Kennedy, and 
A. Kuzmich, Deterministic single photons via conditional quantum evolution, Phys. 
Rev. Lett. 97 (2006), 013601; J. Simon, H. Tanji, J. K. Thompson, and V. Vuletic, In- 
terfacing collective atomic excitations and single photons, Phys. Rev. Lett. 98 (2007), 
183601; C. W. Chou, J. Laurat, H. Deng, K. S. Choi, H. de Riedmatten, D. Felinto, 
and H. J. Kimble, Functional quantum nodes for entanglement distribution over scal- 
able quantum networks, Science 316 (2007), 1316-1320; Z. S. Yuan, Y. A. Chen, B. 
Zhao, S. Chen, J. Schmiedmayer, and J. W. Pan, Experimental demonstration of a 
BDCZ quantum repeater node, Nature 454 (2008), 1098-1 101. 

[36] A. Chia and A. S. Parkins, Entangled-state cycles of atomic collective-spin states, 
Phys. Rev. A 77 (2008), 033810. 

[37] A. S. Parkins, E. Solano, and J. I. Cirac, Unconditional two-mode squeezing of sepa- 
rated atomic ensembles, Phys. Rev. Lett. 96 (2006), 053602. 



24 



Li-hui Sun et al. 



[38] G. X. Li, S. S. Ke, and Z. Ficek, : Generation of pure continuous-variable entangled 
cluster states of four separate atomic ensembles in a ring cavity, Phys. Rev. A 79 
(2009), 033827. 

[39] M. Yukawa, R. Ukai, P. van Loock, and A. Furusawa, Experimental generation of 
four-mode continuous-variable cluster states, Phys. Rev. A 78 (2008), 012301. 

[40] G.-X. Li, Generation of pure multipartite entangled vibrational states for ions trapped 
in a cavity, Phys. Rev. A 74 (2006), 055801. 

[41] R. H. Dicke, Coherence in spontaneous radiation processes, Phys. Rev. 93 (1954), 
99-110. 

[42] Z. Ficek and R. Tanas, Entangled states and collective nonclassical effects in two- 
atom systems, Phys. Rep. 372 (2002), 369-443. 

[43] A. O. Caldeira and A. J. Leggett, Path integral approach to quantum Browian-motion, 
PhysicaA 121 (1983), 587-616. 

[44] C. Cohen-Tannoudj, J. Dupont-Roc, and G. Grynberg, Atom-Photon Interactions, Wi- 
ley, New York, 1992. 

[45] L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, Inseparability criterion for continu- 
ous variable systems, Phys. Rev. Lett. 84 (2000), 2722-2725. 
[46] R. Simon, Peres-Horodecki separability criterion for continuous variable systems, 

Phys. Rev. Lett. 84 (2000), 2726-2729. 
[47] F. Dimer, B. Estienne, A. S. Parkins, and H. J. Carmichael, Proposed realization of 

the Dicke-model quantum phase transition in an optical cavity QED system, Phys. 

Rev.A15 (2007), 013804. 
[48] M. M. Cola, D. Bigerni, and N. Piovella, Recoil-induced subradiance in an ultracold 

atomic gas, Phys. Rev. A 79 (2009), 053622. 
[49] G. X. Li, S. P. Wu, and G. M. Huang, Generation of entanglement and squeezing in 

the system of two ions trapped in a cavity, Phys. Rev. A 71 (2005), 063817. 
[50] D. Kruse, M. Ruder, J. Benhelm, C. von Cube, C. Zimmermann, Ph. W. Courteille, 

Th. Elsasser, B. Nagorny, and A. Hemmerich, Cold atoms in a high-Q ring cavity, 

Phys. Rev. A 67 (2003), 051802. 
[51] B. Nagorny, Th. Elsasser, H. Richter, A. Hemmerich, D. Kruse, C. Zimmermann, and 

Ph. Courteille, Optical lattice in a high-finesse ring resonator, Phys. Rev. A 67 (2003), 

031401. 

[52] J. Klinner, M. Lindholdt, B. Nagorny, and A. Hemmerich, Normal mode splitting and 
mechanical effects of an optical lattice in a ring cavity, Phys. Rev. Lett. 96 (2006), 
023002. 

[53] T. Holstein and H. Primakoff, Field dependence of the intrinsic domain magnetization 

of a ferromagnet, Phys. Rev. 58 (1940), 1098-1 1 13. 
[54] C. M. Caves and B. L. Schumaker, New formalism for 2-photon quantum optics: 1. 

Quadrature phases and squeezed states, Phys. Rev. A 31 (1985), 3068-3092. 
[55] Quantum Squeezing, edited by P. D. Drummond and Z. Ficek, Springer, New York, 

2004. 



